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(N Abstract 



A new detection system KRATTA, Krakow Triple Telescope Array, is presented. This versatile, low threshold, broad energy 
range system has been built to measure the energy, emission angle, and isotopic composition of light charged reaction products. It 
consists of 38 independent modules which can be arranged in an arbitrary configuration. A single module, covering actively about 
t-H 4.5 msr of the solid angle at the optimal distance of 40 cm from the target, consists of three identical, 500 (im thick, large area 
photodiodes, used also for direct detection, and of two Csl(1500 ppm Tl) crystals of 2.5 and 12.5 cm length, respectively. All the 
,i_> ' signals are digitally processed. The lower identification threshold, due to the thickness of the first photodiode, has been reduced to 
0^ about 2.5 MeV for protons (~65 urn of Si equivalent) by applying a pulse shape analysis. The pulse shape analysis allowed also 
I to decompose the complex signals from the middle photodiode into their ionization and scintillation components and to obtain a 
2^ satisfactory isotopic resolution with a single readout channel. The upper energy limit for protons is about 260 MeV. The whole 
, ^ setup is easily portable. It performed very well during the AS Y-EOS experiment, conducted in May 201 1 at GSI. The structure and 
^ performance of the array are described using the results of Au+Au collisions at 400 MeV/nucleon obtained in this experiment. 
O 

■ Keywords: charged particle detection, triple telescope, CsI(Tl) scintillator, large area PIN photodiode, low noise preamplifier, 
^pulse shape analysis 
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1. Introduction 

Charged-particle detection and identification with isotopic 
resolution over a large dynamic range in particle type and en- 
ergy, is mandatory for studies of isotopic effects in heavy-ion 
reactions. Phenomena related to the isotopic composition of 
the reaction system and the emitted products have been shown 
to be useful for exploring the properties of neutron-rich nuclear 
matter as, e.g., encountered in neutron stars [1,2]. Their impor- 
tance will increase with the availability of secondary beams of 
high intensity. Also studies of nuclear structure near the limits 
of stability require increasingly sophisticated and precise detec- 
tion systems. A comprehensive overview of the recent and fu- 
ture developments in this field of instrumentation can be found 
inRef. 0. 

Most of the existing charged particle detectors for the in- 
termediate energy range (up to a few tens or hundreds of 
MeV/nucleon) base their identification on the two- or three- 
fold telescope method EE)- In order to provide the lowest 
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possible identification threshold, the first layer of the telescope 
is usually made of a gas chamber (e.g. DELF [6|, MULTICS 
Q, FASA E), INDRA g), ISiS DID, GARFIELD HQ, FI- 
ASCO Ell) or of a thin Si detector (e.g. FAUST d, INDRA 
13, LASSA El, CHIMERA |I5|. HiRA H3, NIMROD E), 
FAZIA flU). The first, AE, layer is then followed by one or 
two, thicker, Si detectors or scintillators. 

The option with the silicon AE layer, has the advantage of a 
better resolution and is easier to handle, but usually results in 
higher thresholds and is costly. The presented KRATTA mod- 
ules belong to this class of telescopes, however, they have been 
optimized to be budget friendly, without loosing the quality of 
detection. Instead of using the Si detectors of different thick- 
ness, they are using three identical, catalog size, photodiodes 
and two CsI(Tl) crystals. Thanks to the digital signal process- 
ing and the off-line pulse shape analysis, the obtained mass 
resolutions for light charged particles are very satisfactory in 
a broad energy range. On one hand, the pulse shape analysis 
allowed to reduce the energy threshold, resulting from the rel- 
atively thick first layer, by a factor of 3. On the other hand, 
it allowed to effectively double the thickness of the silicon AE 
layer, by combining the ionization components from the first 
two photodiodes, and consequently, to improve the resolution 
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of the first segment of the telescope. 



2. Motivation and requirements 

The main parameters of the KRATTA array have been moti- 
vated by the needs of the ASY-EOS experiment 1 19 1. This ex- 
periment has been designed to study the density dependence of 
the nuclear symmetry energy by measuring flows and isotopic 
compositions of the reaction products from the 197 Au+ 197 Au, 
%Ru+%Ru, and 96 Zr+ 96 Zr react i ons at 400 MeV/nucleon. Dur- 
ing the experiment the most relevant products, neutrons and 
Z=l and 2 particles, have been measured by the LAND [20 1 
detector and the direction and magnitude of the impact vector 
were estimated using the CHIMERA H3 and ALADIN ToF- 
Wall [21] detectors. The KRATTA array has been designed to 
complement the neutron and hydrogen detection with LAND by 
measuring the isotopic composition and flow of light charged 
reaction products up to atomic number Z < 7, and specially, 
to identify the hydrogen and helium isotopes with a resolution 
much better than achievable with LAND. The array was placed 
on the opposite side of the beam with respect to LAND, and 
covered approximately the same solid angle (160 msr). It has 
been designed to detect energetic particles emerging from the 
"mid-rapidity" region of the 400 MeV/nucleon reactions. Mod- 
ular design, portability, low thresholds (below 3 MeV/nucleon) 
and high maximum energy (~260 MeV/nucleon for p and a), 
allow the array to be used in various configurations and experi- 
ments. In particular, it will very well suit the needs of the future 
cyclotron facility at IFJ-PAN in Krakow, planned for proton 
beams from 70 to 250 MeV Last, but not least, the KRATTA 
modules are also compatible with the design of the existing 
Krakow Forward Wall Detector E21 and can be used for reach- 
ing complete coverage of the 2n azimuthal angle. 



3. Active elements and geometry 

The modules of KRATTA are composed of three large area 
HAMAMATSU PIN photodiodes for direct detection J23] and 
of two CsI(Tl) crystals 11241 . The layout and dimensions of these 
active elements are presented in Fig. [T]and their main charac- 
teristics are summarized in Table [T] 



Photodiodes [23 




Csl2 
12.5 cm 



Active Area a 
Thickness " 

Thickness non-uniformity " 

Dead Layers a 

Surface Orientation " 

Full Depletion Voltage b (FD) 

Dark Current @ FD h 

Terminal Capacitance @ FD * 

Rise Time (Laser o pulse) " 



28x28 mm 2 
500±15 urn 
< 3 urn 

1 .5 / 20 um front / rear 

(111) 

120-135 V 

6-16 nA, typ. 9 nA 

190+3 pF 

40 ns 



CsI(Tl) Crystals d 
Tl concentration c 1500 ppm 

Light output non-uniformity e < 7 % 

Shape Truncated pyramids 

Tolerance ±0.1 mm 



Wrapping [26 1 



Reflectance a 
Thickness 



>98% 
65 um 



" Values from technical note. 

b Values from manufacturer's Inspection Sheet. 

c Nominal values. 

Table 1: Main characteristics of the active elements. 



a "Single Chip Telescope", SCT [25 1, configuration and pro- 
vides a composite signal combined of a direct (ionization) com- 
ponent and of a scintillation component coming from the thin 
crystal (Csll). The third photodiode (PD2) reads out the light 
from the thick crystal (CsI2) and, in addition, provides an ion- 
ization signal for particles that punch through the crystal within 
its active area. 



distance from the target [cm] 




Figure 1: Schematic layout of the active elements. 



Figure 2: Single module content. 



The first photodiode (PD0 in Fig. [T| serves as a Si AE detec- 
tor providing the ionization signal alone. It has been "reverse 
mount", i.e. the ohmic side towards the incoming particles. The 
second photodiode (PD1), naturally "reverse mount", works in 



The larger front face of the longer crystal with respect to 
the rear face of the smaller one (see Fig. [T| permits its use 
at an about 2 times larger distance from the target, in config- 
uration that completes the missing half of the FWD phoswich 
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Figure 3: KRATTA in a 7x5 configuration. 



Fragment 




Enf 


F 


>H 


8.3 


89.6 


254.4 


4 He 


8.3 


89.4 


253.9 


7 Li 


9.5 


103.6 


296.5 


20 Ne 


19.9 


231.3 


719.0 


43 Ca 


26.7 


339.7 


1134.2 


91 Zr 


34.0 


513.9 


1911.8 


197 Au 


38.6 


775.8 


3550.9 



Table 2: Lower, E/ olv and intermediate, E,„,, thresholds and up- 
per limits, E up , for selected species (in MeV/u). The thresholds 
correspond to the energy losses in 500 um of Si, in 1000 um of 
Si + 2.5 cm of Csl and in 1000 urn of Si + 15 cm of Csl, re- 
spectively. They have been calculated using the ATIMA tables 
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array [22J. The crystals have been polished and wrapped with 
a highly reflective ESR [26 1 foil, except for the front and back 
windows. The windows have been protected with 6 \im Mylar 
foils. The crystals were optically decoupled. The photodiode 
chips have been glued onto custom-made PCB frames and put 
in close optical contact with the crystal windows. The active el- 
ements have been placed inside aluminum boxes together with 
the charge preamplifiers (see Fig. [2|. The photodiode frames 
and the aluminum housing reduced the geometric acceptance of 
a single module to about 54%. The active solid angle of a mod- 
ule amounts to 4.5 msr. The entrance window has been made 
of a 100 um thick copper foil. During the experiment, 35 mod- 
ules have been arranged in a 7x5 array (Fig. [3J, all placed at 
a radial distance of 40 cm from the target, and operated in air. 
In a spherical coordinate system, with the origin at the target 
and the beam axis defining the reference equatorial plane, the 
columns follow the meridians from the beam level up to about 
27° elevation (latitude) and the rows span the azimuth angles 
(longitude) between 21° and 64° with respect to the beam di- 
rection. The energy thresholds resulting from the thicknesses 



of the consecutive active layers are summarized in Table [2] As 
will be shown later on, the low threshold can be further reduced 
by applying a pulse shape analysis to the ionization signals of 
PD0. The thresholds from Table [2] do not take into account the 
target, the air and the entrance window present in front of the 
active layers. In the case of the ASY-EOS experiment, these 
layers caused increase of the lower threshold, E[ OH ., by about 
12.9 MeV/nucleon for protons and a particles (8.6, 1.8 and 2.5 
MeV/nucleon energy losses due to the 350 um Au target, 40 
cm of air and 100 um of Cu, respectively). For low energy ex- 
periments, the effects of air and the entrance window can be 
minimized by performing the measurement in vacuum and re- 
moving the Cu foil or making it thinner. 

4. Electronics and data acquisition 

The main electronics and data acquisition functions are pre- 
sented schematically in Fig. |4] The photodiodes (3 per module) 
have been reverse biased at 120 V, using the in-house made 
120-channel, remote controlled, high voltage power supplies. 
The signal from each photodiode has been integrated with the 
own-design, low noise, charge preamplifier | 

KRATTA 7x5 mcvp/i 
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Figure 4: Analog, logical and digital signal flow chart. HV PS 
- high voltage power supply. PA - charge sensitive preamplifier. 
V1724 - CAEN digitizers. TRIVA5 - VME Trigger Synchro- 
nization Module l30l . MT - master trigger. FIFO - logical Fan- 
In Fan-Out module. RI04 - VME controller board 1 3 1 1 running 
LynxOS. MBS - Multi Branch System, a GSI acquisition stan- 
dard |[32). SMS - Shared Memory Segment. 

The preamps were supplied with ±6 V and their dynamic 
range spanned about 3.6 V. Three nominal charge gains of the 
preamplifiers have been used, depending on the azimuthal an- 
gle of the module: 44.5, 22.2 and 13.5 mV/MeV, with 1, 2 and 
3.3 pF feedback capacitors, respectively. After optional am- 
plification, the signals have been digitized with the 100 MHz, 
14 bit digitizers [29 1 and stored for the off line analysis. All 
logical and digital electronics modules shown in Fig. [4] have 
been controlled with the RI04 board IT3T1 within a single VME 
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crate. During the experiment, the 14 Flash ADC boards have 
been triggered with an external trigger split and delivered into 
each FADC module. The stored waveforms spanned 5.12 or 
10.24 us (512 or 1024 time bins), with a 2 u.s pre-trigger en- 
abling a precise baseline estimation. The shorter samples have 
been sufficient for the first photodiode supplying the fast ion- 
ization signal alone. The expected data throughput amounted 
to about 5 MB/s, assuming 1 kHz single hit rate. The actual 
data rate did not go beyond this estimate during the experiment. 
The digitizers have been remotely set up and monitored using 
a self-developed software. The data flow has been controlled 
using the standard GSI MBS system Il32l . 



C Qsk \ (RC - tbsXRC - T Fk ) 

(TRS ~ RQ(T RS - T Fk ) (T Fk - T RS )(T Fk - RC) ) 

where, for generalization, Vo is the baseline and At = t - to, 
with to being the beginning of the pulse. 

The two scintillation components |5]l have been introduced 
to account for the fast and slow decay modes of the CsI(Tl) 
crystals. The rise times, trd and Trs, account for the photodi- 
ode, scintillator and the preamp rise times. Overall, the model 
depends on 1 1 parameters listed in Table [3] 



5. Pulse shape analysis 

The pulse shape analysis has two main purposes in the case 
of KRATTA data. First of all, it has to enable the decomposition 
of the signals from the middle photodiode, PD1 (SCT), into the 
ionization and scintillation components. Secondly, it is used 
to resolve masses of particles stopped in the first photodiode, 
PDO, utilizing the relation between the range of a particle and 
the time characteristics of the induced current signal [33 1. 

The following assumptions have been made to accomplish 
these two goals. The preamplifier response has been mod- 
eled using a simple parallel RC circuit approximation |34|: 




v(i) 



iff) _ dVff) Vft) 
~C ~ dt ~RC 



(1) 



where, RC is the feedback coupling time constant, iff) is the in- 
duced current due to the carrier motion in the photodiode, and 
Vft) is the measured voltage pulse. This relation assumes an 
infinite open loop gain, small detector capacitance and a zero 
rise time of the charge integrator, which is an idealization, but 
enables an analytical approach. The induced current has been 
approximated with one direct (ionization), iofi), and two scin- 
tillation, iskft), components, all of the same form: 



i D fi) = Qd 



,-t/TRD 



o-i/tfd 



iskff) = Q. 



Sk 



T RD - T FD 

o-tfrlss _ g-t/TFk 

trs ~ T Fk 



1,2 



(2) 
(3) 



where Q's are the induced charges and the Trd,Trs and the 
t F d, T Fk are the rise and fall times for the respective component. 
The assumed shapes attempt to account for both, the compli- 
cated actual current pulse shape induced by the electrons and 
holes drifting in the photodiode [33], and for the instrumen- 
tal rise time of the preamp. With the assumed preamplifier re- 
sponse ([TJ and the current shapes |2]) and ([3J, it was possible to 
obtain the corresponding analytical model shapes of the wave- 
forms, Vff) = Vq + V D ff) + V s (t), with: 



V D ft) = RC Q D 



-Aj/tjd 



e -M/RC RC 



(RC - t rd )(RC - t fd ) 



TRD 



+ 



(t rd - RC)(t rd - t fd ) (t fd - t rd )(t fd - RC) 



(4) 





Parameter 


Value 




Ionization 






Qd 


Amplitude 




fitted 


Trd 


Rise time 


-90 ns 


fixed/fitted 


T F d 


Fall time 


<300 ns 


fixed/fitted 




Scintillation 






Qsx 


Fast component amplitude 




fitted 


Q S 2 


Slow component amplitude 




fitted 


T F \ 


Fast fall time 


-650 ns 


fitted 


T F 2 


Slow fall time 


-3.2 [is 


fixed 


Trs 


Rise time 


-140 ns 


fixed 




Common 






RC 


Preamp fall time constant 


-6 [is 


fixed 


to 


Time offset 


-2 \is 


fitted 


Vo 


Baseline 




fitted 



Table 3: 11 parameters of the model waveforms. 

The preamp fall time constant parameter RC has been deter- 
mined individually for each chip by selecting the pulses with 
the fast ionization component alone. The resulting RC con- 
stants were smaller than the nominal ones by a few (is due 
to small leakage currents. In order to describe precisely the 
shapes due to particles stopped in the first photodiode, PDO, 
both the time constants, tro and t F o, were fitted and the scin- 
tillation components were obviously not used. In case of PD1 
and PD2, the rise and fall times tro and t ff > were fixed. The 
fits were done using the FUMILI [ 35 1 minimization package, a 
relatively fast and precise implementation of the Gauss-Newton 
algorithm. Constraining some of the parameters was found in- 
evitable, not only to obtain a good description of the waveforms 
fx 2 ) but, at the same time, to maintain the global agreement 
between the reconstructed amplitudes of the three components 
and the predictions of the range-energy tables (see discussion 
of Figs. 
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An additional advantage of using the digitization of the sig- 
nals and the fitting method, was the knowledge of the actual 
charges Q for each component, irrespectively of the substantial 
ballistic deficit (reduction of the amplitude) due to a relatively 
short discharge time of the preamps (RC - 6 (is). For instance, 
the sum Q$ \ + Qsi represents the total light produced in the 



scintillators. For typical values of the parameters presented in 
Table [3] the ballistic deficit amounts to about 5-15% for the 
ionization signal and to about 25 and 50% for the fast and slow 
CsI(Tl) components, respectively. Usually, the fall time con- 
stant of the preamp is a compromise between the level of pile- 
ups, the baseline variation, and the ballistic deficit. However, 
since ballistic deficit is not an issue in our approach, the short 
preamp fall times, of the order of the time span of the waveform 
and of the slow CsI(Tl) decay time, have the additional advan- 
tage of making the maximum and the tail of the pulse visible 
within the digitized sample. 

Waveforms of low energy particles stopped in the first pho- 
todiode, PDO, have been fitted with a single component of the 
form (ffl, with the time constant parameters treated as free fit pa- 
rameters. A more precise time characteristics of the associated 
current pulse (|2| was needed to perform a pulse shape based 
identification of these particles (see discussion of Fig. 13 1. 

The quality of the fits is demonstrated in Fig. [5] The ordinate 
represents the measured voltage, in the FADC channels, which 
is about twice smaller than the actual one due to the matching 
of the 50 Q FADC input impedance. The hits corresponding 
to the selected waveforms are also marked in the identification 
maps presented in the next section. 

Panel A of Fig. [5] shows a waveform for a low energy elec- 
tron, particle or gamma, registered barely above the acquisition 
threshold and stopped in PDO (see also Fig. 13 1. Here the mea- 



sured histogram is visible and the deviations from the smooth 
fit visualize the level of the total noise. The resultant signal 
distortions, including the photodiode, preamplifier, FADC and 
pickup noise sources, have been observed on the level of 0.4 
mV rms, corresponding to about 30 keV. 

Panel B of Fig. [5]presents the quality of the fit with the shape 
given by an ionization component alone Q, applied to parti- 
cles stopped in PDO. Its location in the identification map is 
presented in Fig. 



13 



This fit allowed for derivation of both, the 
rise and fall times and, therefore, also of the mode (position of 
maximum) of the associated current signal 0: 



, trd tfd , TRD 
mode — log 

trd - tfd Tfd 



(6) 



Panel C of Fig. [5] shows the waveform of a high energy a 
particle stopped in the CSI2 crystal. The corresponding hit has 
been marked in Figs. [7|and[T2] 

Panels D-E-F of Fig. [5]show the evolution of the pulse shape 
of an a particle detected and stopped in the SCT as its energy 
increases. The corresponding hits have been marked, when pos- 
sible, in Figs. §[8][TT] 



The fitting of signals from PD2 (Fig. [5] panel C) produces 
small artificial ionization contributions for particles which ac- 
tually do not hit the photodiode (see Fig. [7]). It amounts, on 
average, to 1.7±0.5% of the total amplitude. Correspondingly, 
the artificial contribution of a scintillation component for parti- 
cles stopped in the PD1 photodiode, and thus producing no light 
(see representative hit in panel D), is about 3.8±0.6%. These 
numbers specify the quality of the pulse shape parametrization 
and the systematic uncertainty of the decomposition into dif- 
ferent components. The artificial scintillation component has 
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Figure 5: Example waveforms and their decomposition. Full 
time scale corresponds to 5 (panels A-B) or 10 (panels C- 
F) (is. The ordinate represents the measured voltage in the 
FADC channels. Histogram (when visible): measured wave- 
form. Solid line: sum of all fit components. Dashed line: 
ionization component. Hatched areas: fast and slow CsI(Tl) 
scintillation components. A: low-energy electron, particle or 
gamma ray, barely above the acquisition threshold, stopped in 
PDO (see corresponding hit in Fig. [T3] l. B: a particle stopped 
in PDO (see Fig. 13 1. C: a particle stopped in CsI2 (see Fig. [7] 
T2| ). D: a particle stopped in PD1 (see Figs. 6 8 9|1 1 1 



particle barely punching through PD1 and stopped in Csll (see 
Figs. [6] [HIJTO j JTT|> . F: higher energy a particle stopped in Csll 
(see FigsT[6f]8 [10j[TT]) 



been removed by subtracting its well defined fraction from the 
total amplitude and thus making the ionization and scintillation 
components of the SCT almost perfectly orthogonal. 

6. Performance 



Figures [6]fO] present various identification (ID) maps ob- 
tained by using the parameters of the reconstructed waveform 
components. The strongest lines in Figs. |6p"3"| correspond to p, 
d, t, 3 He, a, and so on, from bottom to top, respectively (see 



Fig. 16 for precise labeling). 

The reconstructed amplitude maps for the first two and the 
last two photodiodes of the KRATTA module are presented in 
Figs.|6]and|7] 

The identification map in Fig. [6] shows a complex spectrum 
with each "ID-line" composed of two parts: an ordinary Si- 
Si hyperbolic segment at low energies, for particles stopped in 
PD1, and a more curved part for particles punching through the 
PD1 photodiode and stopped in the Csll crystal. Due to line 
crossing and its complex structure, this kind of a map cannot 
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Figure 6: AE-E ID map for the first two photodiodes, PD0 vs 
PDl(SCT), for particles stopped in PD1 or in the thin crystal 
(Csll). 
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Figure 8: Decomposition of the map from Fig. [6] PD0 vs scin- 
tillation detected by PD 1 - for particles punching through PD 1 
and stopped in Csll. 
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Figure 7: AE-E ID map of scintillation signals for Csll vs CsI2 
detected by PD1 and PD2. 



be used for identification, however thanks to the presence of 
many characteristic punch-through points and curvatures, it is 
very well suited for energy calibration purposes. 

Fig. [TJshows a AE-E map for scintillation signals detected by 
a second vs third photodiode. Apart from the very good isotopic 
resolution, one can see a substantial background from the sec- 
ondary reactions/scatterings in the crystals and also from parti- 
cles crossing the module at an angle and not originating from 
the target, as well as back-bendings corresponding to particles 
punching through the thick crystal (see discussion at the end of 
this section). Thanks to the pulse shape decomposition, the ion- 
ization component for particles hitting the photodiode (nuclear 
counter effect) at the back of the thick crystal has been easily 
removed. 

Using the individual reconstructed amplitudes for the ioniza- 
tion and for the scintillation components, the map of Fig. [6] 
can be decomposed into the PDO-PDl(Si) and PDO-PDl(CsI) 
components of the SCT segment (Figs. [9] and [8}. This makes 
effectively KRATTA even a four-fold telescope. 
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Figure 9: Decomposition of the map from Fig. [6] PD0 vs ion- 
ization signal from PD1 - for particles stopped in PD1 (produc- 
ing no light). 



The isotopic resolution visible from Fig. [8]can be improved 
by summing up the reconstructed ionization components from 
PD0 and PD1, and thus, increasing the effective thickness of 
the first AE layer to 1 mm of Si. See also inset in Fig. [14] 
Fig. 



1 1 presents a classical AE-E ID map obtained from Fig. 



thanks to the pulse shape decomposition (see discussion of Fig. 
17 for more details on this transformation and Fig. 14 for the 
corresponding mass resolution.) 

Figures[lO]and[T2]show the "Fast-Slow" ID maps for the Csll 
and CsI2 crystals, respectively. The latter represents a variant 
of the standard representation, using the total light vs the ratio 
Slow over Fast (see e.g. |36|). The Fast-Slow representation is 
very useful in many respects in addition to the standard AE-E 
one: One can observe a clear double alpha line (a-a in Figs. [10] 



and 12 1 which is hidden in the AE-E representation behind the 
Li isotope lines. The Fast-Slow map shows a clear "y-line" - 
a strong line below the proton line due to y-rays (see inset in 
Fig. 10 1, which can be precisely isolated and removed from the 
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Figure 10: Fast vs Slow components of the scintillation in Csll. 
The inset shows more clearly a y-line and the p, d, t lines. 
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Figure 11: Summed PD0+PD1 ionization components vs scin- 
tillation in PD1 for particles stopped in Csll. Note a slightly 
improved resolution compared to Fig. |8](see inset in Fig. 14 1 



Figure 12: Total light vs Slow over Fast components of the scin- 
tillation in CsI2. Note, that unlike in Fig. [7] the punch through 
lines do not cross the lower lying isotope lines. 
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Figure 13: Amplitude vs mode of the current signal for particles 
stopped in PD0, obtained from a single Si chip. 



AE-E map. Last, but not least, as can be seen from Fig. 12 the 
punch through segments do not cross the lower lying isotope 
lines in a way they do in case of the AE-E maps (Fig. |7J. Thus, 
the Fast-Slow maps allow to isolate the punch through segments 
in a much more precise way. Unfortunately, the punch through 
segments eventually merge with the y-lines, which makes it dif- 
ficult to discriminate between these two. 



Figure 13 demonstrates yet another powerful feature of the 



pulse shape analysis. It allowed to perform an identification of 
the majority of light particles stopped in the first photodiode by 
plotting the amplitude vs mode of the reconstructed current sig- 
nal (|6}. Due to the constant value of the field within the intrinsic 
region of the PIN diode the enhancement of the resolution for 
stopped particles is not expected to be as pronounced as in the 
case of reverse mount PN detectors [37), nevertheless, the rela- 
tion between the range and the carrier collection time seems to 
be still strong enough to enable the isotopic separation of light 
particles. This method allowed to reduce the lower identifica- 
tion threshold, due to the thickness of the first photodiode, (see 



Table|2]i from 8.3 to about 2.5 MeV for protons, where they are 
still resolvable from deuterons (see two bottom lines in Fig. [T3) >. 
This effectively corresponds to the reduction of the thickness of 
the first AE layer from 500 to about 65 |.im of Si. 



One should stress, that the time scale presented in Fig. 13 
exceeds by about a factor of 3 the collection times estimated on 
the basis of the mobilities of the carriers alone, thus the plasma 
delay, or other effects, seem to be quite important in slowing 
down the collection process. 

The mass resolution for particles stopped in the thin crystal 
(Csll) and in the thick one (CsI2) can be viewed from Figs. 14 
and [15] respectively. 

The background in Fig 



14 resulting mainly from the sec- 



ondary reactions/scattering in the crystal, amounts to about 6% 
for the hits below the p, d, t peaks. For energetic Z = 1 parti- 
cles traversing the thick crystal (12.5 cm of Csl) the secondary 
reaction probability amounts to about 46% (Fig. 15 I. These 



probabilities agree reasonably well with the simple estimates 
obtained using the nuclear collision length of 22.30 cm for Csl 
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Figure 14: Particle identification spectrum (charge distribution) 
obtained from the map of Fig 



11 



for particles stopped in the 
thin crystal. The inset shows the p, d, t peaks (in linear scale) 
obtained from the Fig. [TT](solid) and Fig. [8](dashed histogram), 
respectively. 




Figure 15: Same as Fig. 14 but obtained from the map of Fig. 



|7J for particles stopped in the thick crystal. 



1 38 ], which yields the probabilities of 1 1 % and 43% for 2.5 and 
12.5 cm of Csl, respectively. The background measured under 
the p, d, t lines includes also some contribution from the sec- 
ondary reactions of neutrons and heavier charged particles, as 
well as from the accidental coincidences. The y and punch- 
through hits have been removed from the background, in both 
cases. The high secondary interaction probability obviously de- 
teriorates the identification quality and defines the limits for ap- 
plying the telescope method to intermediate energies. 

7. Discussion and remarks 

Since the parametrization of the pulse shapes is only approx- 
imate and has no deep physical background, one can ask how 
precise is the decomposition into individual components and 
how physical they are. In order to address these questions, the 
AE-E map of the SCT (PDl+CsIl) has been compared to the 
predictions of the ATIMA range-energy tables. Such a compar- 
ison requires the knowledge of the energy calibration of both, 
the silicon photodiodes and of the CsI(Tl) light output. The cal- 
ibration has been performed using the ID map of Fig. [6] (see 
Fig. ~~ 



16 1 which is relatively insensitive to the details of the de- 



composition and is perfectly suited for the energy calibration 
due to its richness in characteristic punch-through points and 
curvatures. 

The calibration routine allowed to adjust the thicknesses of 
the dead and active layers, the energy calibration parameters, 
as well as the light-energy conversion parameters. For the lat- 
ter, a simple integrated Birks' formula (see e.g. |34|) has been 
used, with dE/dx oc A Z 2 /E, yielding a commonly applied two- 
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Figure 16: Same as Fig. [6] but with the superimposed ID lines 
calculated using the ATIMA range-energy tables. 



parameter Light-Energy relation |39|, applicable for particles 
stopped in the scintillator: 



Light = fl] IE - a 2 A Z log 



1 + 



a 2 AZ 2 



(7) 



with E being the energy, A,Z - mass and atomic numbers of the 
stopped particle, 01,02 - the gain and quenching parameters. 

Having this (inverse) calibration it was possible to superpose 
the ATIMA lines on the decomposed ID map for the Single 
Chip Telescope segment of the KRATTA module. The result is 



presented in Fig. 17 and is worth a comment. 




2000 4000 6000 8000 10000 12000 14000 16000 18000 20000 
PD1 Scintillation [channels] 

Figure 17: Decomposed SCT AE-E identification map (ob- 
tained with a single readout channel) with the superimposed 
ID lines calculated using the ATIMA tables. The sequence of 



lines is the same as in Fig. 16 



The observed very good agreement between the recon- 
structed amplitudes and the range-energy table predictions is 
quite non-trivial taking into account the simplicity of the pulse 
shape parametrization. In order to obtain this kind of agree- 
ment, which is still a compromise, it was not enough to set all 
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the parameters free and search for the minimum of the x 1 dis- 
tribution. Such an automatic procedure guarantees the best fit 
to the waveforms, but not necessarily guarantees the follow up 
of the trends defined by the range-energy tables. In general, an 
automatic procedure can lead, depending also on the starting 
values, to some local minima, produce discontinuities or lead 
in directions diverging from the ATIMA lines. The presented 
agreement has been obtained after searching for the optimal 
values and fixing some of the parameters (see Table [3J. The 
agreement with the ATIMA lines could be made even better, for 
instance by freeing the trs parameter, however at the expense 
of loosing completely the mass resolution in Fast-Slow repre- 



sentation (contrary to the case presented in Fig. 10 1. Thus, the 



presented agreement is a result of an iterative procedure of fix- 
ing some of the parameters and also of the energy calibration 
parameters, but once the crucial parameters are constrained the 
fitting proceeds automatically. 

The whole procedure could probably be better automatized 
by using more realistic pulse shape parametrizations, especially 
for individual electron and hole components of the ionization 
signal, taking into account the plasma delay effects, interac- 
tions between carriers, diffusion effects, etc |33|. However, to 
our knowledge, such pulse shape parametrizations are still to 
be developed. Definitely, a more realistic preamplifier response 
function would improve the resolution and would allow for dis- 
entangling between physical and instrumental parameters. A 
more realistic parametrization and analysis would be worth the 
effort of implementing it to e.g. better understand the relation 
between the range and the charge collection times which en- 
able the identification of the particles stopped in the silicon de- 
tectors. Nevertheless, any more sophisticated analysis would 
definitely slow down even more the analysis. 
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Figure 18: AE-E identification map for scintillation signals 
from Csll vs CsI2, with the superimposed ID lines calculated 
using the ATIMA tables. The sequence of lines is the same as in 
Fig. [16] The y-line and punch- through hits have been removed 
(cf. Fig.|7f. 



Finally, Fig. 18 shows the calculated ATIMA lines superim- 



posed on the AE-E identification map for scintillation signals 
from thin vs thick CsI(Tl) crystals. The slight overestimation 



of the AE component for Z = 1 particles results most probably 
from the simplicity of the Light-Energy conversion formula (|7|l 
and (or) from the light output non-uniformity of the crystals. 
The overall agreement is, nevertheless, quite satisfactory and 
the calculated lines can be used not only to derive the energy 
calibration parameters, but also for identification (see Fig. [T5j >, 
after small manual adjustments. 

The energy calibration routine based on the ATIMA range- 
energy tables operated on all three maps (Figs. T6p~8 > simulta- 
neously, allowing for a consistent determination of the calibra- 
tion parameters for all photodiodes and both crystals. Specifi- 
cally, for the CsI(Tl) crystals with 1500 ppm of Tl concentra- 
tion, the quenching parameter «2 was found to be equal to 0.32, 
a value compatible with the one from ll39l and about 20% larger 
than the average one quoted for the INDRA crystals [40]. The 
calibration of the SCT allowed also to estimate the efficiency 
of the scintillation. The combination of a relatively high Tl 
doping, high reflectance of the wrapping and large active area 
of the photodiode, resulted in a relatively high efficiency for 
energy-light conversion. It was found that only about 6 times 
more energy was needed to produce an electron-hole pair in the 
photodiode through a scintillation process in the CsI(Tl) crystal 
than in the direct ionization process in the photodiode. This fact 
is worth noting, since for the "good scintillators" quoted in |34| 
this factor amounts to 15-20. 

Obviously, a drawback of the telescope method at high en- 
ergies, is the high level of the secondary reactions. In order to 
handle this problem, a more sophisticated methods, neural net- 
works and discriminant analysis are being tested, but this goes 
beyond the scope of this article. 



8. Summary 

A new, low threshold, broad energy range, versatile array of 
triple telescopes, KRATTA, has been constructed. The mod- 
ules, equipped with digital electronics chains, allowed for iso- 
topic identification of light charged reaction products. 

Pulse shape analysis allowed for realistic decomposition of 
the complex SCT pulse shapes into individual ionization and 
scintillation components and eventually profit from the, other- 
wise harmful, nuclear counter effect. The isotopic resolution 
obtained using a single readout channel was found to compete 
very well with those obtained using the standard two channel 
readout. The applied pulse shape analysis permitted also the 
identification of particles stopped in the first photodiode and the 
reduction of the identification threshold, due to the thickness of 
the first photodiode, by a factor of three. Thanks to the pulse 
shape analysis, it was also possible to obtain the ballistic deficit 
free amplitudes, which allowed for easy energy calibration and 
identification based on the predictions of the range-energy ta- 
bles. 

The array has met the expectations, fulfilled the design re- 
quirements and performed very well during the AS Y- EOS ex- 
periment at GSI. 
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